Chiral anomaly and r] — rj mixing* 

E.P. Venugopal and Barry R. Holstein 
Department of Physics and Astronomy 
University of Massachusetts 
Amherst, MA 01003 

February 1, 2008 

Abstract 

We determine the rj — rj' mixing angle via a procedure relatively independent of the- 
oretical assumptions by simultaneously fitting rj, rj' reactions involving the anomaly — 
rj,rj' — ► 77, 7r + 7r~7. We extract reasonably precise renormalized values of the octet 
and singlet pseudoscalar decay constants F%, Fq, as well as the mixing angle 9. 
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1 Introduction 



From a strictly theoretical perspective, there exists a significant difference between the octet 
pseudoscalar mesons — tt,K,77 8 — and their singlet counterpart - r?o-0 The former are legiti- 
mate pseudo-Goldstone bosons whose masses vanish in the chiral limit while the latter is not 
due to the anomalous breaking of the axial U(3) symmetry down to SU(3). However, in the 
real world, this does not seem to make much difference. Indeed, the physical eigenstates — 
r], rf — are mixtures of octet and singlet components 

77 = 7/8 cos 9 — T)o sin 9 

rj' = rjs sin 9 + t)q cos 9 (1) 

and the mixing angle 9 is an important quantity in confronting theoretical calculations with 
experimental results in these systems. 0] 

The mixing angle can be evaluated in various ways but a standard procedure involves 
the diagonalization of the 77, rj' mass matrix, which at lowest order yields a mixing angle 
9 ~ —10°. This can easily be seen by writing the mass matrix as 

m 2 = ( m i m | ) (2) 

Employing the Gell-Mann-Okubo (GMO) relation to fix m\ as0 

m l = \{^ 2 K -ml) (3) 

and diagonalizing, one can determine 771q 8 , m\ and 9, in terms of m^m^, yielding 

9 = -9.4°, m 2 8 = -0Mm 2 K , m Vo = 0.95 GeV (4) 

However, the GMO relation, Eq.[| is valid only at lowest order — 0(p 2 ) — in chiral per- 
turbation theory. The inclusion of 0(p 4 ) corrections to the relation results in significant 
changes. Characterizing these via 

ml= l -{Am 2 K -ml){l + 5) (5) 



a leading logarithm estimate [|j] 



4 2 

O m K !„ m K 
Z (47rJV) 2 111 (471-FV) 2 

4m 2 K — m\ 



(6) 



yields 5 ~ 0.16, which yields 

9 £ -20°, m 2 8 = -O.8I774, m = 0.90GeV (7) 

suggesting a doubling of the mixing angle. A full one-loop chiral perturbation theory calcu- 
lation confirms this finding. It is also interesting that this solution is consistent with the 
assumptions of simple U(3) symmetry wherein ?7 8 , 770 have the same wavefunction, leading to 

9 2v / 2 /771 — 771 s \ r, . . 

m 2 = — ~ -0.9 m\ 8 

08 3 \m + m s ) K v ; 
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Note that this result is strongly dependent upon chiral symmetry breaking effects on the 
GMO prediction. 

An alternative and independent approach to the problem involves the use of the chiral 
anomaly, which is responsible for the well-known 71,77,77' — > 77 decays. In the case of 
7T° — ► 77, at lowest order — 0(p 4 ) — the anomalous (Wess-Zumino-Witten) chiral lagrangian 
predicts 0] 



A~o. 



»77 



3ttF 



2/3 



(9) 



where £j, hi are the polarization, momenta of the outgoing photons, and F is the pion decay 
constant in the chiral limit. A leading log calculation of the chiral corrections reveals that the 
dominant effect is simply to replace F by its physical value F n = 92.4MeV.|§ The resulting 
amplitude is guaranteed by general theorems to remain unchanged in higher chiral orders. M 
One then finds that the predicted amplitude 



Fjryy (0) 



3ttF w 



0.025 GeV~ l 



is in excellent agreement with the corresponding experimental value |T0] 

F wyy (0) = (0.024 ±0.001) Gel/" 1 



(10) 



thus providing the confidence that one may be able to analyze the T], rj' decays with a similar 
precision. 

In case of the rj, rf — > 77 decays, one must include both, mixing as well as renormalization 
of the octet, singlet couplings, yielding the predicted amplitudes Q 



-^7/77(0) 



8 Fn / 



aN c /F 

3y/3irF n \F 8 

37S;U sm ' + 2 ^ cos ' 



(12) 



From the PDG one extracts experimental values PD 



FfjjjiO) 

Frfy ; (0) 



0.0249 ±0.0010 GeV- 1 
0.0328 ± 0.0024 GeV' 1 



(13) 



In order to solve the system, however, we require an additional assumption since there 
are three unknowns — Fg, Fq, 9 — and only two pieces of data. The usual approach in this case 
is to use the leading log prediction of one- loop chiral perturbation theory to predict B 



m 



K 



In 



m 



K 



mz 



(4^)2 (4^)2 (AnF n 



1.30 



(14) 



1 Note that we implicitly assume here, as do other workers, that all 0(m^, m^L/A^.) effects, where A x ~ 
A-kFk is the chiral scale, are included in the renormalization of the pseudoscalar couplings Fg, Fq and in the 
mixing angle 9. This does not have to be the case, but appears to be borne out by the consistency of the 
results obtained from treatments of differing manifestations of the anomaly, as we show below and as others 
have found. 
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and then solve for F Q ,6, yielding 

^ = 1.04, 9 = -20° (15) 

It is interesting to note that these results are quite consistent with those obtained from 
the one-loop analysis of the mass matrix — i.e. 9 « —20° and F /F n consistent with the 
value — unity — which one would have if the singlet state and the pion were to have the same 
wavef unction. 

While this agreement is satisfying, the extraction of these mixing parameters requires 
certain theoretical inputs, either Eq. [14| or Eq. |6], and it is interesting to inquire whether 
one can predict the mixing angle purely phenomenologically. As we shall show, the answer 
is yes, provided one utilizes the additional information available from the anomalous decays 

V,V' -> 7T+7T-7.Q 

In the next section then we show how the decays rj, rf — > 7r + 7r _ 7 can be analyzed in order 
to isolate the chiral anomaly This involves a careful study of final state interactions and 
unitarity constraints in order to realistically extrapolate to zero four-momenta, as required 
by the anomaly. In the concluding section, we apply these results to evaluate 9, F s , F in an 
essentially model independent fashion. 



2 Analysis of 77, rf — > tt + tt 7 Decays 

Our goal is to use the experimental data on r], 7/ — > 7r + 7r~7 in order to isolate the value of 
the anomaly in these decays. The resulting numbers can then be fit to the appropriate the- 
oretical expressions, thus allowing extraction of the renormalized mixing angle and coupling 
constants. To this end, we define 

A ri: n'^iTiT'Y = Br),T)'( s i s ttit) ^ VOL ^ tykvV+aV-fl (16) 

where p±, k are the outgoing 4-momenta, e M is the photon polarization vector, s = rni„, and 
s n7T = (p + +p-) 2 . The chiral anomaly (cf. Figure 1(a)) requires 

B -' (0 '°^l5fe(t sme + ^f CO8( ') (17) 

Note, however, that the chiral anomaly, strictly speaking, only constrains the form factors 
B^^'is, s^) at zero four momentum — 5^ ^/(0,0) — while the experimental input occurs at 
s = m ^,V' s tttt — 4m^. One indication of this fact is that, using the simple energy- independent 
form given in Eq. [H] to calculate the decay rate, one obtains for the ^-channel the value 
r^-TrTpy = 35.7 eV compared to the experimental rate of r^™, = 64 ± 6 eV. For the 
channel things are even worse — r.^/_ vr7r7 = 61 ± 5 KeV while the theoretical value arising 
from use of the simple anomaly amplitude is a factor of twenty less! We conclude that 
in order to extract values for the anomaly in these transitions, it is absolutely essential to 
correctly model the energy dependence of the amplitude in the physical region. 
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rj/rf 




7 



(b) (c) 
Figure 1: Shown are contact (a) and VMD (b,c) contributions to 77,77' — > 7r + 7r~7 decay. 



This problem is not unique, of course, to the 77, 77' — > tttt^ system and has been addressed 
previously in the extraction of the anomaly in 7 — > 3n from the Primakoff-effect data of 
Antipov et al.|12| In this case, one also has a clear prediction of the chiral anomaly. |13 
Writing 



.4 



7— >3-7r 



F 37T (s,t,u) e 



fiua/3 



where s = (p + + p~) 2 ,t = (p + + po) 2 , u = (p_ + po) 2 , the chiral stricture demands 

eN c 



^(0,0,0) 



12ttFJ 



9.7 GeV 



(18) 



(19) 



The effects of p-wave pi-pi interactions at low and moderate energies are known to be reason- 
ably described by vector dominance. [0] In particular, the models described in ref. |[L5| , which 
incorporate vector dominance and chiral symmetry, when applied to the 7 — > 37r reaction, 
provide a form 



F 37T (s,t,u) 



1 



^(0,0,0) 



m. 



m. 



t ml 



u 



(20) 
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which matches on to the anomaly at zero four-momentum but also offers a plausible extension 
into the physical region. This form must be modified, of course, in order to confront real data 
since unitarity demands the presence of branch cuts and consequent imaginary components 
in the form factor. This is clear from a one-loop chiral perturbation theory calculation, which 
yields pi 



37T 



where 



s,t,u) 



F 37r (0,0,0) 



1 3m 2 
1 H + 

2m 2 p 



m; 



24vr 2 F 2 



3 77T 

— In — + F(s) + F(t) + F(u) 

4 mi 



(21) 



F(a) 



2(1 



4ml 



s 

4m? 



-4m| j n 1+ 



-1+ 



4m?— s , _1 

— 21 — tan 



- 2 s > 4m! 



(22) 



4m?- 



2 s < Ami 



Here we note that the imaginary component of the function F(s) is given in terms of the 
energy-dependent width of the rho meson via 



ml 



24vr 2 F 2 



ImF(s) 



m, 



(23) 



where 



r P (a) 



9pTTTT S 

A8irm n 



'I _ 4m M3/2 



(24) 



This one-loop form is no doubt appropriate in the near threshold region. However, once 
Snw > 10m 2 or so, one does not anticipate that a simple one-loop description will be adequate. 



In order to address these difficulties, the use of an N/D form has been suggested. ||17|| In this 
approach, one utilizes the Omnes function, which encodes information concerning the pi-pi 
interaction |18| 

' s f°° ds' 5i(s') \ 



DAs) 



exp —- 



71 



It 



(25) 



where Si(s) is the (1 = 1) p-wave pi-pi phase shift at center of mass energy ^/s. There are 
two ways to proceed at this point. 



i) One can use the experimental phase shifts, with some assumptions made about their 
asyptotic form. In our case we took the values quoted by Froggatt and Peterson, JT9 
which are given up to ^Js = 1 GeV and assumed a constant value after that. We label 
the Omnes function obtained in this way as Dl xp (s). 



ii) One can employ a simple analytic form[| 

Dx(s) = 



96tt 2 F2 



fc4 

mi 



m; 



24tt 2 F 2 



F{s) 



(26) 



which has been shown to provide an approximate description of the empirical tttt p- 



wave phase shifts in the low energy region. [51] We label the Omnes function obtained 
via this procedure as Df na,l (s). 
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Using either of these forms, and postulating an N/D form of the 7 — > 3tt amplitude as 

1 m? m? m? 

F 3w (s,t,u) = --,l, ;: (0);i ' " " " " 




m 2 p — t m 2 p — u ' 



x f£- fz. (27) 



one can see that Eq. ^ matches onto the one-loop chiral result in the limit of low energies 
and onto the vector dominance form when unitarity inspired logarithms are dropped. 

Whether such an N/D form accurately describes the data for 7 — > 3tt awaits the arrival 
of sufficiently precise information from CEBAF and CERN. However, it certainly appears to 
satisfy the various criteria which nature demands, and suggests the treatment of the related 
rj, rf — > 7r + 7r~7 decay amplitudes in a parallel fashion. 

In the case of the rj, rf — ► 7r + 7r~7 decays we can proceed similarly. In this case, the 
one-loop chiral perturbation theory calculation gives 

1 1 m 2 4 20 3 

Bl- l °° P (s, s™) = B v (0, 0)[1 + ((-4ml + 3^) In ^ + ^F(s^) - -ml + — s^) 

P (28) 

while vector dominance (cf. Figure l(b,c))yields 

^(3,3 CT ) = ^(0,0)[1 + | *™ ] (29) 

Certainly, in order to treat the decay of the //, one must go further and include unitarity 
effects via final state interactions. One very simple approach is to include the (energy- 
dependent) width of the rho-meson in the propagator via 

(30) 



m p s n7T Tnp Smr im,pV p(s- K - K ) 

This use of vector width-modified vector dominance already makes an important difference 
from the simple anomaly — tree level — results (especially in the case of the rj') , changing the 
predicted decay widths from the values 35 eV and 3 KeV quoted above to the much more 
realistic numbers 

62.3eV, T 

rf — 7T7T7 

Ql.bKeV (31) 

if the parameters 

F 8 /F w = 1.3, F /F V = 1M, 6 = -20° (32) 

are employed. However, this approach does not reproduce the one-loop chiral form in the 
low energy limit. 

In order to determine a form for the final state interactions which matches onto both the 
one-loop chiral correction and to the vector dominance result in the appropriate limits, we 
postulate an N/D structure, as in the related 7 — > 3tt case — 



Byj — 7T7T7 (S, S7T7r) -B?; — 7T7T7 (0, 0) 



1 + asi, 
1 — c + c- 



F>i(s 7Tn ) 



(33) 
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where for the Omnes function we use one of the two forms itemized above and a, c are free 
parameters to be determined. In order to reproduce the coefficient of the F(s 7r7r ) function, 
which contains the rho width, we require c = 1. On the other hand, matching onto the VMD 
result at 0(p 6 ) can be achieved by the choice a = l/2m 2 p . Thus in the case of the 77 the form 
is completely determined. Since the 7/ spectrum is closely related and is dominated by the 
presence of the rho we shall assume an identical form for the rf case. Using these forms we 
can then calculate the decay widths assuming the theoretical values for the anomaly. Using 
the parameters given in Eq. |32 one finds, for example, 



i)Dt 
ii)D\ 



cxp / 



anal / 



77 — 7T7T7 



r/ — 7T7T7 



65.7eV, 
69.7eV, 



1 T\' — 7T7T7 
1 r/'— 7T7T7 



m.2KeV 



(34) 



There is a tendency then for the numbers obtained via the analytic form of the Omnes 
function to be somewhat too high. 




100 150 
E 7 (MeV) 

(a) 




100 150 
E y (MeV) 

(b) 



Figure 2: Shown is the photon spectrum in rj — > 7r + 7r~7 from Gormley et al.|2(| as well as 
various theoretical fits. In the first figure, the dashed line represents the (width-modified) 
VMD model. The (hardly visible) dotted line and the solid line represent the final state 
interaction ansatz Eq.[3^ with use of the analytic and experimental version of the Omnes 
function respectively. The second figure shows the experimental Omnes function result 
(solid line) compared with the one- loop result (dotdash line). 



We can also compare the predicted spectra with the corresponding experimentally de- 
termined values. As shown in Figure 2, we observe that the experimental spectra are well 
fit in the 77 case in terms of both the N/D or the VMD forms, but that the one-loop chiral 
expression does not provide an adequate representation of the data.[^] In the case of the 
corresponding 77' decay the results are shown in Figure 3, wherein we observe that either 
the unitarized VMD or the use of N/Dl xp provides a reasonable fit to the data (we get 
X 2 /dof=32/17 and 20/17, respectively), while the use of the analytic form for the Omnes 
function yields a predicted spectrum (x 2 /dof=104/17) which is slightly too low on the high 
energy end. However, for both 77 and 7/ we see that our simple ansatz — Eqf33] — provides a 
very satisfactory representation of the decay spectrum. 
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400 600 800 

m™ (MeV) 

Figure 3: Shown is the photon spectrum in r/ — > 7r + 7r~7 from Abele et al.|[24|| as well as 
various theoretical fits. As in Figure 2, the dashed line represents the (width-modified) VMD 
model. The dotted and solid lines represent the final state interaction ansatz Eq.^ with 
use of the analytic and experimental version of the Omnes function, respectively. Here the 
curves have been normalized to the same number of events. 

3 Evaluation of 77 — rf mixing parameters 

Our conclusion in the last section was that if the mixing angle and pseudoscalar coupling 
constants were given values consistent with present theoretical and experimental leanings, 
then the predicted widths and spectra of both 77, 77 — > 7r + 7r~7 are basically consistent with 
experimental values. Our goal in this section is to go the other way, however. That is, 
using the assumed N/D forms for the decay amplitude, and treating the pseudoscalar decay 
constants F 8 , F as well as the 77 — 77' mixing angle 9 as free parameters, we wish to inquire 
as to how well they can be constrained purely from the experimental data on 77, 77' — > 77 and 
77, 77' — > 7r + 7T - 7 decays, with reasonable assumptions made about the final state interaction 
effects in these two channels. 

On theoretical grounds, one is somewhat more confident about the extraction of the 
threshold amplitude in the case of the lower energy 77 — > 7r + 7r~7 system. Indeed, in this case 
the physical region extends only slightly into the tail of the rho unlike the related 77' decay 
wherein the spectrum extends completely over the resonance so that there exists considerable 
sensitivity to details of the shape. Thus a first approach might be to utilize only the two- 
photon decays together with the 77 — > 7r + 7r~7 width in order to determine the three desired 
parameters. In this fashion one finds the results shown in Table 1. We observe that the 
results are in agreement, both with each other and with the chiral symmetry expectations — 
F 8 /F n ~ 1.3, F /F n ~ 1, and 9 ~ —20°. However, the uncertainties obtained in this way are 
uncomfortably high. 
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F /F v 


e 


VMD 

N/D anal 
N/Df p 


1.28 ±0.24 
1.49 ±0.29 
1.37 ±0.26 


1.07 ±0.48 
1.02 ±0.42 
1.02 ±0.45 


-20.3° ±9.0° 
-22.6° ±9.6° 
-21.2° ±9.3° 



Table 1: Values of the renormalized pseudoscalar coupling constants and the 77 — 77' mixing 
angle using the 77, — 77 and 77 — 7r7r7 amplitudes in a three parameter fit. 





F B /F V 


F /F„ 


9 


VMD 

N/D anal 
N/Df p 


1.28 ±0.20 
1.48 ±0.24 
1.38 ±0.22 


1.07 ±0.04 
1.09 ±0.03 
1.06 ±0.03 


-20.8° ±3.2° 
-24.0° ±3.0° 
-22.0° ±3.3° 



Table 2: Values of the renormalized pseudoscalar coupling constants and of the 77 — 77' mixing 
angle obtained from a maximum likelihood analysis using the 77, 77' — 77 and 77, 77' — nn^f 
amplitudes. 



In order to ameliorate this problem, we have also done a maximum likelihood fit including 
the 7/ — 7T7T7 decay rate, yielding the results shown in Table 2. We observe that the central 
values stay fixed but that the error bars are somewhat reduced. The conclusions are the 
same, however — substantial renormalization for F 8 ~ l-3F n , almost none for F ~ F n , and a 
mixing angle 9 ~ —20°. These numbers appear nearly invariant, regardless of the approach. 



4 Conclusion 

Before summarizing the results of our above analysis, it should certainly be emphasized that 
we are not the first to undertake the program of isolating the anomaly from the 77, 77' — 7m r f 
data. Indeed, there has been considerable work in this regard, both on the theoretical 



side |TT| , [23f as well as experimentally, including the most recently published 7/ — > 7r + 7r~7 
data. [24] The recent analysis of ref.|23|] leads to results quite different from ours in both the 
mixing angle as well as the renormalization of the pseudoscalar couplings. On the other 
hand, ref. 0] (at least for the model labeled Mi) finds a somewhat smaller mixing angle 
(9 ~ —16°) and pseudoscalar renormalization (F 8 /F ~ 1.1). 

However, there is an important difference between these analyses and our own. In refs. p3[ 



and |23| , the decay amplitude is written in terms of a piece due to the anomaly (parametrized 
by Ex, X = 77, 77') and a component due to the rho pole (parametrized by Fx, X = 77, 77'). The 
ratio Ex I Fx is then fitted, through a minimization procedure, to produce the experimental 
spectrum and partial widths. In our analysis, the parameters of the two pieces are fixed 
a priori to reproduce the results of one-loop chiral perturbation theory [|l(| (by fixing c = 1) 
and VMD |15|] (by fixing a = l/2m 2 p ). (Indeed a recent analysis of other 1=1 7r + 7r~ processes 
found that only within a model such as ref.[l5|], which links chiral symmetry and VMD, 
could the data be fit consistently |25|]). However, following the picture of ref.[|15j, we do not 



10 



include a non- resonant coupling in 77r + 7r (as considered in models Mi and M 3 of ||23|| ). Our 
successful fits of the experimental data speak for themselves. 

In previous treatments of the 77, 7/ system via the anomaly, which have omitted rj, i]' — tttt^ 
constraints, the mixing angle 9 has generally been determined only at the cost of theoretical 
assumptions about the renormalization of the octet pseudoscalar coupling constant F 8 with 
respect to F n . We have in this paper asked whether it is possible to obtain the mixing angle 
in a fashion relatively independent of such theoretical assumptions by simultaneously fitting 
77, rf — > 77 as well as 77, 77' — > 7r + 7r~7 decays. As shown above, the answer is affirmative. 
However, one must incorporate some sort of model for the final state interactions of the pi-pi 
system in order to extrapolate down to zero four-momentum where the anomaly obtains. 
We have argued that the N/Dl xp form given in Eq. |33| is reasonable both on theoretical 
grounds — matching both the requirements of VMD and of low energy chiral symmetry — and 
via successful fitting of the experimental spectra. Using this form and the PDG values for 
77, 77' — > 77 and rj, 77' — > nnj amplitudes we have obtained values 

F 8 /F w = 1.38 ±0.22 F /F n = 1.06 ±0.03 9 = -22.0° ± 3.3° (35) 

which are quite consistent with those obtained in previous analyses which required assump- 
tions about chiral symmetry breaking. One can then assess these results in two different 
ways. Although it is our contention that the assumptions made above concerning pion-pion 
interactions are relatively model-independent and that the numbers given thereby in Eq. |35| 
are quite solid, one could also take a contrary view that the forms utilized for final state 
interactions do require critical dynamical assumptions. In this case, however, we would argue 
that via three quite different routes 

i) mass matrix analysis including GMO breaking 

ii) 77,77' — > 77 analysis with assumptions made about F 8 /F n 

iii) simultaneous analysis of 77, 77' — > 77 and 77, 77' — > 7r + 7r~7 with (minimal) assumptions 
concerning final state pi-pi interactions 

one finds virtually the same value of the mixing angle — 9 ~ —20° — and for pseudoscalar 
couplings — F 8 /F n ~ 1.3, F /F n ~ 1.0. In any case, we would assert that these values are 
now strongly (and independently) confirmed from within the chiral anomaly sector. 
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